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We present a more detailed picture of the infrared regime of Landau gauge Yang-Mills theory. This 
is done within a novel framework that allows one to take into account the influence of finite scales 
within an infrared power counting analysis. We find that there are two qualitatively different infrared 
fixed points of the full system of Dyson-Schwinger equations. The first extends the known scaling 
solution, where the ghost dynamics is dominant and gluon propagation is strongly suppressed. It 
features in addition to the strong divergences of gluonic vertex functions in the previously considered 
uniform scaling limit, when all external momenta tend to zero, also weaker kinematic divergences, 
when only some of the external momenta vanish. The second solution represents the recently 
proposed decoupling scenario where the gluons become massive and the ghosts remain bare. In this 
case we find that none of the vertex functions is enhanced, so that the infrared dynamics is entirely 
suppressed. Our analysis also provides a strict argument why the Landau gauge gluon dressing 
function cannot be infrared divergent. 

I. INTRODUCTION 

The combined effort of functional approaches and lattice gauge theory led in the past years to a comprehensive 
picture of the qualitative features of the infrared (IR) limit of Yang-Mills theory in Landau gauge |T}{3]- This qualitative 
information is encoded in a set of IR power laws for the Green functions of the theory. These can incorporate important 
aspects of the confinement mechanism for gluons within the scenarios of Kugo-Ojima ^4^ and Gribov-Zwanziger |S] 
and serve as a basis for the inclusion of matter fields and the challenging problems of quark confinement [BHH] and 
spontaneous chiral symmetry breaking ^ . 

Previous studies of the Dyson-Schwinger equations (DSEs) [TJ [SJ [TU] suggested that the qualitative aspects of the 
Yang-Mills IR fixed point structure in Landau gauge are already known |llfll4| . It had been shown that when 
the DSEs and functional renormalization group equations are combined, within the class of scaling fixed points, 
where no masses are dynamically generated, there exists a unique solution |15j . This conclusion was possible by the 
complementary constraints obtained from the two different hierarchies of equations and we will show here that, taking 
into account the constraints provided by the existence of the skeleton expansion, this is already implied by the DSE 
system alone. Yet, analog to the propagators, which depend on a single external scale, it was expected that the IR 
limit of a Green function is uniquely determined by a single IR scaling law. In the case of the vertices, however, 
the situation is more diverse. Besides the appearance of a multitude of different tensor structures, the corresponding 
form factors are also functions of several distinct momenta. Therefore, there are in general different combinations 
of momenta that can become soft and lead to IR divergences. In present studies the implicit assumption was that 
vertices become IR divergent if and only if all scales go to zero. This led to general results for the IR scaling of 
arbitrary vertices in this uniform limit jTl] that were later extended to arbitrary dimension jTB]. Here we show that 
this picture - although qualitatively correct - needs to be refined in the sense that there are additional kinematic 
singularities that characterize the IR-regime. 

It is known since the early work of Taylor [T7] that the IR limit of the ghost-gluon vertex in Landau gauge is not 
infiuenced by radiative corrections when an external ghost momentum vanishes. This result is used as the starting 
point for various IR analyses and is explicitly confirmed by our analysis. In contrast we find a kinematic singularity 
in the 3-gluon vertex. Although this kinematic singularity which refiects the sensitivity to ultrasoft gluon exchange is 
rather mild, ~ (p^) with k > 0.5 |T21|T3], it could be conceptually important. First, due to to the parametrically 
larger support in loop integrals kinematic singularities should have a sizable impact on the quantitative results for 
Yang-Mills Green functions. Further, it has been argued recently that the Slavnov- Taylor identity for the 3-gluon 
vertex suggests that the gluon propagator is IR divergent [18]. This argument relied on the assumption that the 3- 
gluon vertex is finite when only a single momentum vanishes. Our results show that this is not the case for the scaling 
solution and thereby the corresponding conclusion cannot be drawn. In addition we show here that the DSEs for 
the gluonic vertices lead without any approximation or assumption to the condition that the gluon dressing function 
cannot be IR divergent jTS] as argued by Mandelstam [12]. Finally and even more important, the techniques developed 
here allow one to study the quark sector of QCD, where a corresponding kinematic singularity in the quark-gluon 
vertex is induced [6] . In contrast to the gauge vertex where the kinematic singularity is induced semi-perturbatively, 
cf. |20| . the one in the quark-gluon vertex is much stronger divergent due to a self-consistent enhancement mechanism. 
Thereby it induces long ranged gluonic interactions that can confine quarks [51 [7] . 

In contrast to the assumptions of present scaling analyses the DSE system entails also the possibility of dynamical 
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Figure 1: DSE for the ghost-propagator. Thin and thick lines and dots represent bare and proper propagators and vertices. 
AU proper Green functions are amputated and the external lines are shown only for illustration purposes. 




Figure 2: DSE for the gluon-propagator. 



mass generation. In addition to the IR scaling solution a decoupling solution where the gluons become massive 
has recently been suggested pT1 - E5] . see also E5] . This solution is also observed in recent lattice simulations 
[3l • Whereas there is only a single scaling solution the decoupling case seems to allow a whole family of solutions 
characterized by the IR gluon mass that continuously connects with the scaling solution. These different solutions are 
obtained depending on the chosen boundary conditions for the functional equations [23_, cf. also [21]. In case of the 
decoupling solution certain IR Green functions are not dominated by IR modes but by finite scales of the order of the 
characteristic scale Aqcd- We extend the power counting formalism in order to take this possibility into account in 
our IR analysis. Thereby we confirm that the decoupling is compatible with the vertex equations and indeed presents 
a second solution of the full DSE system. However, we find that within the decoupling solution none of the vertices 
is IR enhanced. Correspondingly, no long ranged interaction is reflected in the Green functions of the fundamental 
local degrees of freedom. 

This article is closely related to another article PU] on the IR limit of Yang-Mills theory and the behavior of the 
vertices. There we provide explicit, analytic solutions for the 3-point vertices in a semi-perturbative approximation 
and confirm and elaborate on the results discussed here. 



II. GENERALIZED INFRARED FIXED POINT ANALYSIS 

In this section we present the general formalism that allows one to study the non-perturbative IR scaling behavior 
of Green functions. In particular, we take into account that the IR fixed points can be influenced by large scales, 
like external momenta of Green functions that remain flnite in the IR limit and can lead to kinematic singularities, 
or dynamically generated masses. Our analysis is based on the system of Dyson-Schwinger equations that describe 
the complete non-perturbative dynamics of the considered theory. The equations for arbitrary vertex functions can 
be derived algorithmically as shown in [55] and as has been implemented in the Mathematica package DoDSE. The 
algorithm can be represented by diagrammatic replacement rules applied to the underlying equations for the 1-point 
functions and is sketched in Appendix [K\ 




Figure 3: First version of the DSE for the ghost-gluon vertex. The first order in a skeleton expansion consists of the first three 
diagrams on the right-hand side. The fourth diagram, which features the IPX ghost-gluon scattering kernel, contributes only 
at higher order. 
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Figure 4: Second version of the DSE for the ghost-gluon vertex. Here the leading order in the skeleton expansion of the right 
hand side consists of the first three diagrams in the first and second lines, respectively. 




Figure 5: Full DSE for the 3-gluon vertex. To leading order in a skeleton expansion the last graphs in the first respectively 
third line are absent. The restriction to the f-loop graphs reduces to the leading order in ots when dressed Green functions are 
replaced by their tree-level expressions. 

The equations for the propagators and primitively divergent vertices of Landau gauge Yang-Mills theory for instance 
are given in figs. [l]to[6] As can be seen, the DSEs for the primitively divergent Green functions are in general not 
closed. Instead they are part of an infinite hierarchy of coupled equations and thereby it might seem hopeless to 
make any definite statements about existence and uniqueness of possible fixed point solutions. However, as has been 
demonstrated in fl4], this is even possible on the level of a mere IR power counting analysis. The working hypothesis of 
such an analysis is that in the IR regime the theory can still be described by the field content of the classical action. It 
involves in the present case as a tool the concept of a skeleton expansion of higher order correlation functions in terms 
of skeleton diagrams. This means that all vertex functions that do not appear in the Lagrangian are replaced by a loop 
expansion that features only primitively divergent vertices and propagators. Consequently the system of equations 
is closed, as it consists only of the DSEs of the primitively divergent vertex functions. The difference compared to a 
perturbative expansion of these higher vertex functions is, however, that all the primitively divergent Green functions 
in the expansion loops are dressed. Since the DSEs employed here are based on one-particle irreducible (IPI) vertices, 
graphs that involve vertices that are not primitively divergent do not contribute to leading order in the skeleton 
expansion. E. g. for the vertices of Landau gauge Yang-Mills theory the leading order in the skeleton expansion 
is discussed in the captions of figs. |3] to [6] The skeleton expansion yields a closed DSE system for the primitively 
divergent Green functions but involves graphs of arbitrary loop order. In the framework of a power counting analysis 
it is possible to assess the IR scaling of all these graphs since higher orders can be generated by simple extension 
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Figure 6: DSE for the 4-gluon vertex - for conciseness we only show the 1-loop contribution in the skeleton expansion, 
reduces to the leading order in as when dressed Green functions are replaced by their tree-level expressions. 
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rules [TJ. As discussed in more detail below, such extensions by additional loops should in particular not increase the 
degree of IR divergences since higher order diagrams would otherwise be more and more divergent and the description 
in terms of the underlying degrees of freedom would break down. Correspondingly the lowest order in the expansion 
already entails in general the leading IR scaling. The solutions obtained with a skeleton expansion are by construction 
solutions of the full system of DSEs, where the scaling of higher order Green functions is already determined by the 
leading order skeleton graphs. Note that the skeleton expansion does not miss additional solutions where the non- 
trivial behavior is triggered by n-point functions, n > 4 |7]. In the case of scaling fixed points the IR analysis has 
already been performed without the tool of the skeleton expansion in ^5, 27 . Furthermore, the equivalence between 
the constraints obtained from the skeleton expansion and constraints derived from functional renormalization group 
equations has been shown in jUj. Hence the use of the skeleton expansion is well justified. 

Whether a solution, where the underlying degrees of freedom are still valid in the IR regime, actually exists can only 
be decided by an explicit analysis of the non-perturbative dynamics. Once found, a corresponding solution exists 
then independently of the skeleton expansion and directly verifies the working hypothesis. Yet, our analysis does not 
exclude other possible solutions where this is not the case. A consequence of our working hypothesis is that there 
should be no singularities in the external momenta in two- and three-point functions as long as all scales are finite. 
Such a singularity arises e.g. in the perturbative running of the strong coupling in the form of the Landau pole, but 
in this case it merely refiects the insufficient dynamical treatment. Yet, in a fully non-perturbative analysis it would 
signal the breakdown of the description in terms of the fundamental gauge degrees of freedom. In contrast, it is for 
instance generally expected that the physical degrees of freedom in Yang-Mills theory are glueball bound states which 
are reflected as singularities in higher gluonic n-point functions with n > A. These singularities arise at finite invariant 
momenta corresponding to the masses of the glueball states. In Euclidean space these poles are on the negative half 
axis of the squared momentum variable and do not interfere with our analysis. 

The motivation for an IR analysis is that certain theories, including in particular the important case of gauge theories, 
are scale invariant at the classical level. Yet, the quantum theory can dynamically generate a scale A via dimensional 
transmutation. Whereas e.g. in a perturbative analysis of Yang-Mills theory the corresponding scale Aqcd is 
renormalization scheme dependent but otherwise uniquely defined by the Landau pole, non-perturbatively there is 
no unique definition and it can be chosen as an arbitrary hadronic scale. For the scaling solution of Landau gauge 
Yang-Mills theory [llj, discussed in more detail below, it could for instance be identified with the scale where the 
gluon propagator reaches its maximum. The leading behavior of all Green functions far below the scale A should 
by renormalization group arguments be described by a power law scaling with appropriate IR exponents. Although 
such a power law form is obtained analytically in conformal field theories and has been observed in scale free regimes 
of simpler models, for Yang-Mills theory the existence of an IR power law scaling has to our knowledge not been 
rigorously proven, yet, and presents the main assumption we make in our analysis. To be precise, in this work we 
determine the IR solutions of Yang-Mills theory that feature IR power law scaling but cannot exclude other possible 
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solutions that involve an IR momentum dependence with a more complicated analytic structure. In particular it is 
conceivable that in addition to a power law scaling there are additional logarithmic dependencies, like it is the case in 
the UV regime of the theory. We do not explicitly take into account this possibility, but we note that such additional 
logarithmic dependencies would only refine the scaling and not change the dominant power law structure. Thereby 
the qualitative results of this analysis remain unaltered. 

Whereas the propagators depend only on a single external momentum scale and correspondingly have a unique IR 
behavior, the vertices involve several momenta and can feature different IR power laws in different kinematic sections. 
To expose the IR behavior, we perform a tensor decomposition of a given proper Green function Ty with appropriate 
Lorentz indices fii and internal indices Oj in some tensor basis Tt that is analytic in the momenta 



i^v)l'^.'.'.ll^{qi,--- ,9„) = ^r.„,t (g?,gi-g2,-- • ,ql) (TtT^'^.'.'.^l^ (qi, ■ ■ ■ (1) 

t 

so that the corresponding dressing functions T^^t depend only on scalar arguments. Under the assumption of IR 
scaling. Green functions can in principle have distinct IR power laws whenever any combination of the external 
momenta vanishes. The dressing functions generally take this power law scaling form only in the deep IR regime 
p| <C A. For intermediate momenta of the order of A they can involve a more intricate structure and in the UV 
region they feature a different asymptotic scaling behavior, that is e.g. in Yang-Mills theory given by the perturbative 
logarithmic RG running. Nevertheless, the full dressing functions can be formally decomposed into a sum of terms 
that include the individual IR power laws in the different kinematic regimes 




r r 2 2\ 1/2 2\ I t'l V^l' ' \ ln\ 

r«,t(9i,gi-g2,--- ,9„) =2^7«,t(9i,9i-'72,--- ,g„) ( — ^ — p -] ■ (2) 

In the following we explain this expression in detail. Here, the set of functions pf defines the scaling variables in 
the different kinematic limits I which are chosen analytic in all arguments and identical for all Green functions. The 
Sl ( are the corresponding power law exponents. In the IR limit that the scaling variable vanishes, the corresponding 
function 7^ j in the decomposition eq. (|2j) is defined such that it depends only on finite momentum ratios 

7^,,t(9?, 91-92, ••• ^^7^,t Q^l^ • (3) 

where the scaling variable pi depends only on the momenta qi. that vanish in the kinematic limit. In particular if 
the scaling variable depends only on a single external momentum, 7^ j becomes constant in the IR limit. In general, 
the IR exponents 6^ ^ also depend on the specific tensor. In our analysis we will not distinguish between the different 
tensor structures and are only interested in the IR exponents of the most singular dressing functions = mint t) 
which will dominate in the IR. 

We should stress, that the additive decomposition eq. ^ which explicitly separates the different IR divergences of 
a dressing function into different summands does not represent an IR approximation and is exact on all momentum 
scales. Yet, it is surely far from unique. This ambiguity is reflected in the functions 7^ ^ which include the proper mid 
and high momentum behavior and also depend on the choice of the arbitrary scale A. This means that the explicit 
IR power law is compensated in the functions 7^ ^ outside of the IR region. As discussed above, for n-point functions 
with n > 4 they may involve singularities when all combinations of external scales are finite but they are by definition 
regular whenever external scales vanish. The only requirement on the partition is the complete decomposition of 
possible IR singularities. For the determination of the IR power law exponents, however, the detailed form of the 
decomposition at finite scales is irrelevant as discussed below. 

An important IR limit is determined by a uniform scaling variable that is given by a function p„ that fulfills 



Plill--- :<ll) ^0 ^ gi , • • • , g„ ^ A qf/pl,--- , ql/pl constant . (4) 
Such a function is for instance provided by the "Euclidean norm" 
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Figure 7: The unique IR singular kinematic configuration of tfie propagator (left) and the different IR singular kinematic 
configurations of the 3-point functions (right). 
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Figure 8: IR singular kinematic configurations of the 4-point functions. In this case there can even be kinematic singularities 
when all external momenta are hard. 



Since in this scaling limit all external momenta scale uniformly to zero this presents the fully scale-invariant case. 
Previous studies assumed that all Green functions are finite in other kinematic limits and thereby have only considered 



this "conformal"^ case [TH - flB] . We will discuss this important approximation in section IV 

In general, however, there are various different kinematic limits where the pi can also depend on any subset of the g;. 
As we will argue below in section |V] in Yang-Mills theory such a more general IR behavior is realized that involves 
kinematic singularities according to eq. ([2| when only a subset of the external momenta vanishes. We will denote those 
external momenta that scale to zero as soft and those that stay fixed when the IR limit is taken as hard momenta. 
In particular, all soft external momenta are assumed to be much smaller than A in order to ensure the applicability 
of IR power laws. By momentum conservation it is impossible that exactly n ~ 1 momenta of an n-point function 
become soft, but e.g. all other cases when 1 < i < n — 2 momenta tend to zero can define a separate scaling limit. 
The simplest such kinematic limit is when a single external momentum becomes soft and all others remain hard where 
the corresponding scaling variable in eq. ^ is pf = qf. For each distinct external field this presents a different limit 
described by a separate power law exponent. These present the only kinematic limits in case of 3-point functions, 
as illustrated in fig. |7] As shown in fig. [8] for the 4-point functions there are already two distinct limits when one 
or two external legs are soft. In addition, there can even be IR singularities for 4-point functions when all external 
momenta are hard but differences of momenta become small corresponding to the exchange of a soft momentum in 
the respective intermediate channel. For higher order Green functions the number of IR exponents rises further. 

Now let us discuss how to determine the IR behavior of the integrals arising from the loop corrections in the DSEs. 
Inserting the parameterization eq. ^ into a general 1-loop integral Q yields a sum of different integrals of the form 

S = HE A"^^"^'"^'^"^^ liPh^'^M ({(^ + n + Q^f '-''pl {{ik + Qf} , {Q'}Y' (6) 



where m and v run over the internal propagators and vertices in the loop and ly over the the different IR sensitive 
limits of the vertex v. The Qi are linear combinations of the external momenta and the functions 7^" as well as any 
other analytic terms arising from the tensor decomposition are combined in the kernels They may involve 

additional power law divergences in the external momenta alone but are analytic in the loop momentum. Such 
integrals arise in addition to the explicitly shown IR sensitive kinematical cases also for the individual terms from the 
tensor decomposition eq. (jlj. We note in particular that due to the scaling form eq. ^ the dependence on the scale 



The term conformal refers in our context only to the scale invariance property but not the additional symmetries included in the 
conformal group. 
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A factors out of the actual integral. 

The scalar integrals that remain after an appropriate tensor decomposition are in general far too complicated to be 
performed analytically. However, in the special case that the vertices are constant explicit solutions are known. The 
1-loop two-point integrals are given by the simple analytic form 

{2nf^ ' ^' -I ^ ' r(-i.i)r(-z.2)r(d + i.i + j/2) ^ ^' 

for d/2-|-i^i + 1^2 < 0, d/2 + i^i > and d/2 + 1^2 > 0, where the expression is convergent. An analytic expression for 
the corresponding IR 3-point integrals in terms of hypergeometric functions is known PSI and we discuss how to 
extend this result to the Euclidean regime of these integrals in [2CP . 

However, in order to solely determine the IR scaling of a given integral an explicit solution is not required. In the 
limit that some external scales tend to zero and others stay finite one has a clear scale separation. As is well known 
from effective field theories and shown in detail in appendix |B] this allows one to effectively decouple these scales. 
The basic steps are pictorially illustrated in fig. |9] An integral involving both hard and soft scales is first divided 
at an arbitrary intermediate scale. The now limited support of the two integrals allows one to Taylor-expand the 
integrand so that in the IR limit, where the lowest order is exact, the individual integrals depend only on either soft 
or hard scales and the respective other scales factor out of the integral. Next, the two integrals are extended over 
the whole scale range again. This leaves two correction terms that can be expanded once more and combined to a 
scale-independent integral extending over all scales. It merely presents a counterterm to cancel possible divergences 
introduced in the process of separating the initial integral. Thereby the initial integral is effectively decomposed into 
two separate integrals that depend only on one of the two classes of scales respectively, whereas in the soft part the 
dependence on the hard scales factors out of the integral. On dimensional grounds such integrals must scale as a 
function of this single scale and the power law exponents can then be determined by a mere power counting analysis. 
This works since the integrals are effectively dominated by the poles of the integrand that are in turn determined by 
the external momenta. In particular, these distinct integrals can each feature a different IR scaling behavior. Thereby, 
when assessing the IR scaling of a given correlation function via its DSE, all these partial integrals over the different 
kinematic regions of the initial loop integral can be analyzed individually and the most IR divergent term determines 
the IR scaling. As discussed in appendix |B] it can be necessary to repeat this decomposition when multiple soft and 
hard scales are present. After a complete decomposition the different kinematic contributions of the initial integral 
can be displayed and treated just as distinct Feynman graphs, as shown in fig. [TO] for the soft-particle limit of the 
3-point function, where the momenta of the internal propagators in the individual contributions are restricted to be 
soft (s) respectively hard {h) as indicated. 

Furthermore, it is necessary to consider the presence of mass scales in the IR analysis. In addition to the case of 
theories that have a sector which involves explicit mass scales at the classical level which is coupled to a conformal 
sector, like QCD, there is even the possibility of dynamical mass generation in theories that are scale-free at the 
classical level, like Yang-Mills theory studied below. In general a massive behavior of a (scalar) propagator D is 
defined by 



lim D (p ) = const . (8) 



This yields for the corresponding IR exponent 6 
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In this case the IR scaling is broken since an explicit mass scale 

= lim (p^) (10) 
is present and the propagator can be parameterized e.g. by the alternative forms 
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Figure 9: Pictorial representation of the decomposition of loop integrals described in detail in appendix[B] A loop integral that 
involves both a soft s and a hard scale h can in the limit s <ti h he divided at an arbitrary intermediate scale and decomposed 
into two independent loop integrals and a counterterm that cancels possible divergences introduced in this process. Each of 
these two integrals involves only a single scale whereas the other scale factors out of the corresponding integral. This allows a 
direct power counting of the IR scaling behavior of the individual contributions. 




Figure 10: IR singular kinematic regions in the decomposition of the triangle integrals contributing to the 3-point functions 
in the soft particle limit. There are two different kinematic regions that involve soft loop momenta corresponding to the two 
inequivalent ways to route the hard momentum through the loop. 



In this form it is clear that integrals can also be dominated by scales of the order of the mass, as discussed in detail 
in appendix |B] In the presence of mass scales there are relevant contributions from corresponding hard momentum 
scales in the loop integration even in the case that all external scales are soft, as shown in fig. [TT] If these hard 
contributions dominate such a dynamically generated mass presents a self-consistent solution of the DSE system and 
poses an additional possibility for the IR behavior that has to be taken into account. 

In the uniform limit the scaling is very different in the massive and massless case. When there are only massive 
propagators in a loop integral the integral is entirely dominated by hard scales even in the uniform limit and does not 
scale with the external soft scales, see appendix [B] When only some propagators are massive the integral can again 
be decomposed and receives contributions from both soft and hard momentum regions, as is again shown in appendix 
|B]and is depicted in fig. [TT] In contrast, in the massless case the integrals are IR dominated in the uniform limit 
and the contribution from hard loop momenta is strongly suppressed. This can be seen by expanding the analytic 
functions /C{/„} in eq. (joj) as well as the arising propagators in the limit qi <C k, where k is the loop momentum. 
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Figure 11: Possibly dominant kinematic regions of 1-loop integrals of the propagator (left) and the 3-point function in the 
uniform limit (right). The labels s and h denote soft momenta that vanish when the scaling variable tends to zero and hard 
momenta that stay finite in the limit, respectively. 
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The leading order term is indeed independent of the external momenta. In case there are no other scales involved, i.e. 
in particular a 7^ 1, cf. appendix [B] this scale independent integral as well as the third term is canceled or removed 
in the renormalization process, as discussed below. The linear term vanishes likewise in the symmetric integration so 
that only the last term remains and the contribution from hard modes is actually suppressed in qf/k'^- Yet, when 
there is any explicit scale in the integral, like a hard external momentum or a mass (corresponding to a = 1), the 
leading term in the expansion eq. ( 12 ) is not canceled and there is no suppression. In order to treat the different 



cases within a common framework and allow for all possible solutions, including those where masses are dynamically 
generated, we count in all cases the contributions from both the hard and soft momentum region of the loop integral. 
In case a mass is generated dynamically it arises from quantum fluctuations in the "mid-momentum regime", i.e. hard 
scales of the order of the generated mass. To establish whether a massive solution actually exists and to obtain a 
quantitative value for the mass requires a numerical solution of the functional equations for all momenta. However, a 
mere IR analysis can nevertheless determine if a dynamically generated mass presents a possible dynamical solution 
of the DSE system. This is the case since the above expansion eq. (12 1 shows that the different cases of a massless 



and a massive propagator are effectively discriminated by their IR exponents. In the power counting study below we 
implement this qualitatively different scaling of the propagators arising in loop integrals via the symbol /i, where the 
index i stands for the corresponding particle species and which is defined by 



I for 5i = 1 , i.e. a massive IR behavior 

fay S,^l (^2) 



Since eq. ( 12 ) shows that the contribution from the hard momentum regime of a loop integral is suppressed only 
if all propagators are massless, the suppression of the whole loop integral is then given by adding the product of 
the Hi of all propagators in the loop to the loop's IR exponent. We stress that although in the considered region of 
hard momenta the dressing functions do not feature the IR scaling form, the above prescription is meaningful for the 
power counting analysis since the massive behavior is nevertheless reflected in the anomalous IR power laws by the 
IR exponent Si ~ 1. 

The loop integrals in the DSEs are UV divergent and must be renormalized. Since the leading Green functions have 
positive mass dimension power law divergences seem to pose a threat to a naive power counting analysis because 
hard momenta could contribute even in the case that all external momenta are soft. Actually, the above suppression 
explicitly shows the known fact that all power law divergences, like e.g. the quadratic divergences in the propagators, 
are momentum independent in scale invariant theories. If present, these identically cancel between the different loop 
graphs in a DSE. Whereas this is generally guaranteed by the perturbative UV behavior so that the contributions from 
asymptotically large momenta /c 3> A to the integral cancel, in the IR scaling limit the cancellation of contributions 
from hard modes is according to eq. ( |12[ ) already realized for loop momenta k ^ p. Similarly there are dedicated 
counterterms that by construction explicitly cancel any additional singularities introduced in the above decomposi- 
tion, cf. appendix iBl This leaves at most logarithmic momentum-dependent divergences to be removed within the 
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renormalization procedure. As shown e.g. in [TTl [T^ this is possible multiplicatively by corresponding wave-function 
and vertex renormalization factors that multiply the corresponding tree level terms as well as the loop corrections 
in the DSEs. Since the logarithmic momentum-dependence is not considered within the power counting analysis 
anyhow, it is clear that the problem of renormalization does not interfere with the IR scaling analysis as could have 
been expected considering that they operate in opposite limits of the loop integrals. In principle the analysis can also 
be performed directly for already subtracted DSEs where the issue is omitted from the outset. Thereby, in general 
UV divergences and the detailed renormalization procedure do not affect the IR analysis. Yet, there is one important 
exception to this since under certain conditions discussed below the renormalization prescription can determine the 
boundary conditions for the DSE system [22 . 

After a complete decomposition of the loop integrals in the DSEs, the IR scaling of each contribution is determined 
via power counting of the anomalous and canonical dimensions of the involved Green functions. The most divergent 
contribution - corresponding to the minimal IR exponent - dominates and determines the scaling of the Green function 
on the left hand side of the DSE. Finally it is important to note that there can be various cancelations caused by the 
symmetries of the considered theory that require a careful analysis. This will be discussed in detail in section [V| for 
the case of Landau gauge Yang-Mills theory. 

III. DISPROVAL OF THE NAIVE INFRARED SLAVERY SCENARIO 

After this general section we will now begin our IR analysis of Yang-Mills theory. The equations for the propagators 
of Landau gauge Yang-Mills are given in figs. [T]and[2] These equations have been studied extensively with appropriate 
ansaetze for the vertices. The corresponding DSEs for the primitively divergent vertex functions of Yang-Mills theory 
are given in figs. |3]to|6] For the ghost-gluon vertex there are two qualitatively different versions derived via the ghost- 
respectively gluon-part of the path integral. The leading order in a skeleton expansion is given by the two triangle 
diagrams that have been analyzed in |30] within a semi-perturbative analysis, where it was found that this vertex is 
hardly changed from its tree-level form. The other vertices have so far been discussed only via IR scaling analyses 
[T4l - fl6| which we will detail in this work. Explicit IR results for the 3-point vertices are given in [20J. 
Before we study the implications of the system of equations for the IR limit of Yang-Mills theory, we will demonstrate 
the power counting method explicitly for an important example. Consider the first two diagrams of the second line 
in the DSE for the 3-gluon vertex fig. |5] The contribution from these two graphs to the right hand side of the DSE is 

g.) = i / ^ {iTo)Z% [k - q,) {k^q„ -k, q,) Df^ (fc) 

+ (ro)';t/^ Df^ {k - q{)V^^\l (k - gi, -fc, gi) B% (fc)) , (14) 

where D represents the dressed gluon propagator, F the dressed vertices distinguished by their indices, and Fq the 
bare versions. We note that each of the two integrals depends on only one of the two independent external momenta, 
because the bare 4-gluon vertex is momentum independent. Therefore, the two integrals cannot exactly cancel each 
other for general momenta. With the uniform IR scaling exponents 6gi for gluon propagators and the corresponding 
uniform exponent for the 3-gluon vertex each of the integrals as well as its sum scales in the uniform limit as 

AFf;,(gi,<Z.) ((P^)"^'")' (Z)^^'^"' - (p2)H^-3.+2.., _ ^^^^ 

We see that the canonical scaling of the 3-gluon vertex given by the \ drops out since it appears both on the left 
and right hand side of the DSE and as expected it is sufficient to count only anomalous IR exponents in the uniform 
scaling limit. Each of the contributions on the right hand side of the DSE - or several of them - could dominate and 
determine the scaling of the 3-gluon vertex on the left hand side of the DSE. The leading term is the one with the 
smallest IR exponent. Correspondingly, the IR exponent of the two exemplary graphs has to be larger or equal to the 
left hand side in case they dominate. This leads to the very general condition 

6lg<8lg^2&gi ^ 6gi>0. (16) 

The only other possibility would be that the two integrals are canceled identically by other diagrams in the DSE 
which correspondingly would have to have precisely the same kinematic dependence. We pointed out before that the 
two graphs have a very special kinematic structure given by 
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(91-92) = F^t", (gi) + iq,) . (17) 

Evidently all other graphs involve the two momenta in a manifestly non-linear way - e.g. already the propagators in 
the loop induce manifest non-linearities - and thereby cannot have the above simple property. For instance consider 
the case of the third gluon loop graph that includes the proper 4-gluon vertex. Its integral representation is 

, f d'^k d f , , , Z{k + q,) {k + qi)^{k + qi)p\ i.^ ^ W /^A 

' J (^ ^ ~(M^ ii^' ) -'^^ ^+^^1' ^ 

(18) 

where the gluon dressing functions Z generally involve non-integer powers. Even in the simplest case dgi — 1 where 
the propagators become trivial the above property would be in contrast to the IPI nature of the proper 4-point 
vertex. Since the momentum dependence of the other contributions in the DSE is even more non-linear it is fair to 
conclude that there are no identical cancelations between the individual diagrams in fig. |5] This yields the direct 
constraint that the gluon dressing function cannot be singular as found in the Mandelstam approximation |19] . The 
same relation is obtained from the gluon-loop corrections given by the last three diagrams in the second line of the 
4-gluon vertex DSE fig. |6] Note that this result is a direct prediction of the full vertex DSEs and does not involve 
any assumption or approximation. Moreover, it is independent of the detailed renormalization prescription used for 
the DSEs. Interestingly, all studies that found an IR enhanced gluon propagator made uncontrolled assumptions on 
the form of the 3-gluon vertex, whereas such a non-enhanced IR behavior of the propagator was indeed found as soon 
as the vertex equations were considered dynamically within simplified semi-perturbative analyses, cf. |10) . 



IV. CONFORMAL SCALING 



In this section we will perform an IR fixed point analysis for Landau gauge Yang-Mills theory but will first neglect 
the possibility of kinematic singularities or dynamically generated mass scales. Thereby we consider the idealized, 
conformal case that all external momenta as well as the considered correlation functions scale with a single uniform 
scaling variable p„, as defined in eq. (|4|, and are finite in other kinematic limits. In this case eq. (pi reduces to 




r„ (91,91-92, ••■ ,g„j ~7i, ('7i,9i-'72,--- ,9„j ■ 72 • (19) 



Correspondingly, we neglect the possibility that contributions from hard loop momenta are relevant in the IR limit 
and assume that only the IR regime dominates all loop integrals. As will be shown in the next section this assumption 
is too simplified and the solution structure obtained in this section is incomplete. 

Since the Dyson-Schwinger equations form an infinitely coupled system of equations it is necessary to reduce it to 
a manageable form, as discussed in sec. [IT] This is done via a skeleton expansion that presents a loop expansion 
in terms of dressed primitively divergent correlation functions. The skeleton expansion can be generated from the 



leading order graphs by a finite set of extensions that increase their loop order [14] as shown in fig. 12 As a necessary 
condition for the skeleton expansion these extensions must not increase the IR exponents for a given Green function 
since otherwise successive extensions would make it arbitrary singular. Thereby, assuming that the skeleton expansion 
is not explicitly divergent, provides additional constraints for the IR exponents of the primitively divergent vertices. 
When there is only a single uniform scaling variable p„, so that all external momenta scale with it, qf ^ up 
to possible subleading corrections an integral has to scale as a power of it (or possibly as a logarithm in case the 
naive power counting yields a constant). As noted before, in four dimensions all canonical momentum dependence 
cancels and it turns out that it suffices to count anomalous powers of to assess the IR behavior of a general loop 
correction. The leading dynamical contribution on the right hand side of its DSE determines the scaling of a given 
Green function. With the corresponding uniform IR exponents for the ghost propagator Sg^, the gluon propagator 
5gi, the ghost-gluon vertex 5gg, the 3-gluon vertex 5^g and the 4-gluon vertex S^g we can analyze the IR scaling limit 
of the DSE system for the five primitively divergent Green functions. Here we skip the superscript u on the uniform 
exponents as it appears in eq. (19) for better readability in this section. From figs. [T]and|2]one reads off the power 



counting relations for the IR exponents of the propagators^ 



^ The tadpole contribution in the gluon DSE does not appear here, since it does not depend on the external momentum. 
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Figure 12: Possible extensions of given graphs to generate higher orders in the skeleton expansions. The crossed out propagators 
are part of the initial graph and are not counted. 



-6gh = (0, 6gg + Sgli + Sgl) , 

-Sgi = min (0, S^g + 2dgi,Sgg + 25gh, 2(533 + 4<5gi, (549 + iSgi) 



and correspondingly for the vertex functions from figs. |3]to[6] 



(20) 



&gg = min (0, 2Sgg + 25 gh + Sgl,53g + Sgg + Sg,,. + 25 gl) , (2l) 

= l^in (0, 25 gg + 25 gh + Sgl,25gg + 5gh + 25gl,35gg + 25 gh + 3(5g/, (^Sg + 25gg + 5 gh + 4:5 gl) , 

= min(0, 25gg + 35gh, 25zg + 3(5g;, (^Sg + 25gu5Ag + 25gl,'i5:ig + ^5gu5Ag + (53g + 4:5 gl) , 

(54g = (0, 35gg + 4:5gh,353g + 4:5gl,5ig + 25gl,253g + 35gl,5ig + ^Sg + 35gl,4:53g + 65gl,54g + 253g + b5gu25ig + 45 gl) 

Besides the constraint 5gi > obtained before there are additional analogous constraints from the linear terms in 



the vertex equations (21 1. These constraints are weaker than those obtained from the corresponding RG equations 



studied in [151 127] and are not sufficient to ensure a unique solution of the system of equations. Therefore, we will 
not exploit them in the following and thereby circumvent the problem of possible cancelations discussed above. 



In its present form the nonlinear system of the DSE conditions (20 1 and (21 1 is rather complicated due to the involved 



minimum functions. However, if one requires that there should exist a sta 
in fig. [12] yield the additional much stronger constraints 



ale skeleton expansion, the extension graphs 



25gg + 25gh + 5 gl > Q , 
^3g + 5gg + 5gh + 25gl > , 

2<53g + 35gi > , 

5ig + 25gl > . 



(22) 



The combination of IR exponents in the first two of these constraints is precisely the one arising from the two triangle 
diagrams in the first equation for the exponent of the ghost-gluon vertex eq. (21 1, so that this equation becomes 
unique whereas the second one remains non-trivial 



A 5a 



min (0, 25gg + 25 



gl 



^gh 



(23) 



Since both of these equations have to hold simultaneously it is clear that the second equation has the trivial solution, 
too. To see that it is not important which DSE we start with, consider only the second equation and assume that 



the second term dominates. In this case the ghost-gluon vertex would be given by 5. 



gg 



-25, 



gl 



Ogh- 



Plugging this 
that 5gi is 



III 



in the first constraint from the skeleton expansion, eq. (22 1, yields 5gi < 0. Since we know from section 
non-negative, only the trivial solution 5gi = 5gh remains. 
This non-renormalization condition has previously been used as a starting point in the analysis |il4) . It was supported 
by the gluon transversality in Landau gauge but required the additional assumption that the ghost-gluon scattering 
kernel should not be strongly divergent. Here it arises directly from the physical requirement of a stable skeleton 
expansion. The two other constraints in (22 1 precisely remove the non-linearities in the equations (21 1 for the gluon 



vertices which could lead to a self-consistent enhancement of these equations. Inserting this IR scale independence of 
the ghost-gluon vertex into the other equations and using the constraints eqs. (22 1 as well as the previously shown 
relation 5gi > yields 
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Sgh = min (0, Sgh + Sgi) , 
-5gi = min(0, 2(5g,,) , 
(53g = min (0, i5gh) , 

5^g = min (0, Ugh^iS^g + 4,5g/) . (24) 

As it stands the naive system of equations above has only the trivial solution that all anomalous IR exponents vanish 
since for 5gi > the first equation yields directly 6gh = and the rest follows trivially. 

However, so far it has not been taken into account in our scaling analysis that the DSEs have to be renormalized. As 
shown in |11||TB] when a propagator is divergent, it is possible to do this renormalization at = which cancels the 
tree level part identically. Thereby, the renormalization condition can act as a boundary condition for the DSEs, as 
observed in [23]. Since the exponent of the gluon propagator is non-negative, as shown before, this renormalization 
prescription is not possible there. In contrast due to this positivity the mixed loop correction in the ghost equation 
proves that the ghost propagator can indeed be divergent and the two exponents are then connected by 6gh = —Sgi/2. 
Because of this connection the equation for the gluon propagator is dominated by the ghost loop, but it becomes 
trivial and does not determine Sgi. The solution depends therefore on a free parameter 6gi = 2k > 0, 6gh = — t < 0. 
With the expressions for the propagators this immediately gives S^g = — 3k from the ghost loop which shows that the 
ghost loop dominates the gluon loop also in the four gluon vertex DSE and we have d/^g — —4k. Therefore, there is 
only a single scaling solution of this system depending on a real parameter k > when the relations arising from the 



condition of a stable skeleton expansion eq. (22 1 are taken into account. This is the solution obtained previously in 
[T3HTB] where it was shown that the above ghost dominance mechanism holds for arbitrary n-point functions. Ghost 
loop contributions to gluonic correlation functions and minimal mixed loop contributions to correlation functions 
involving ghosts dominate and all gluonic corrections are suppressed due to the scaling of the propagators despite the 
strongly divergent gluonic vertex functions. Here, in particular the trivial solution 6i — , Vi, obtained for a generic 
renormalization prescription is also contained in the above scaling solution for k = 0. The restriction to solutions 
of the scaling type arises here since we implicitly assumed that the integrals are always dominated by modes in the 
vicinity of soft external momentum scales. In the next section we will see that this is too simplified and that the 
result can be changed by the presence of large scales. 



V. INCLUSION OF SOFT SINGULARITIES AND MASSES 



As discussed in detail in section [IT] the IR counting is complicated by the possibility of kinematic divergences or 
dynamical mass generation. In these cases there are hard external scales present that do not tend to zero when taking 
the IR limit and care has to be taken to assess how a certain correlation function scales with the soft momenta. 
Whereas for the propagators there are unique anomalous IR exponents, in general there can be different IR exponents 
that describe how a correlation function scales in different kinematic sections. The uniform limit where all external 
scales go to zero uniformly presents the conformal case studied above. Beyond this, there are for the 3-point functions 
possibly distinct IR exponents when only a single momentum vanishes, as illustrated in fig. [7j According to eq. Q, 
we discriminate these different kinematic limits by upper indices and denote the corresponding uniform indices in the 
following by d^g and S^g. There are in addition the exponents for the ghost-gluon and 3-gluon vertices, when a gluon 

momentum vanishes, 5^g and 5^g, respectively, and for the ghost-gluon vertex, when a ghost momentum vanishes, 
(59^. As illustrated in fig. [s] for the 4-gluon vertex, studied in the uniform limit in [31_, there can be two distinct 

additional kinematic exponents S^g when one or S^g when two gluon momenta vanish. In addition there could even 
be IR divergences when all external momenta are large but differences of momenta become small. Because of the 
symmetry of the vertex this can be described by one additional IR exponent 5\g corresponding to the exchange of a 
soft momentum in an intermediate channel. 

Although possible, it is rather cumbersome to determine the IR exponent of the 2-loop graphs appearing in some DSEs 
taking into account the various possible kinematic singularities. Fortunately, from what we know from the analysis in 
the uniform limit this should not be necessary. There the observed strong ghost dominance strongly suppressed gluonic 
contributions compared to the leading ghost loops. Because of the absence of primitively divergent 4-point interactions 
involving ghosts all of these leading contributions involve only 3-point vertices and are 1-loop graphs. Assuming that 
the additional kinematic singularities do not entirely change this property motivates a truncation scheme involving 
only dressed 3-point vertices neglecting all 2-loop graphs and those involving dressed 4-point vertices. The only graphs 
that involve bare 4-point vertices that are not 2-loop are those discussed in sec. |III| and which provided the mere 
constraint Sgi > 0. Therefore we will not have to discuss them here again. The truncation we study here is also similar 
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to the one obtained from a 3-loop expansion of a 3PI action analyzed in [3? but with the difference that in the DSEs 
there is one bare vertex in every graph. 

To recapitulate, the general assumptions of the IR analysis introduced in section [II] are 

(i) that the local degrees of freedom are valid to describe the system, 

(ii) that the IR regime of the DSE system can be analyzed using a skeleton expansion and 

(iii) that the IR is described by power law scaling. 
Furthermore, in this section only 

(iv) we truncate the system to 1-loop graphs to simplify the analysis. 

We emphasize that once solutions are found it can be explicitly verified that they present consistent IR solutions 
of the full DSE system. Green functions are then dominated by the leading order skeleton contributions involving 
primitively divergent Green functions, i.e. the assumptions (i) & (iv) are indeed fulfilled. The assumptions (ii) & (iii) 
merely limit the class of solutions that can be found with our technique. 

In the case that several external scales are present the loop integral can receive relevant contributions from fiuctuations 
in the vicinity of all these different scales. Therefore, as discussed in section [TTj it is necessary to decompose the 
momentum integral into distinct integrals which in the IR limit are entirely dominated in the respective soft or hard 
momentum regions. These integrals involve only one characteristic scale each and their IR scaling can then again 
be determined by power counting. Furthermore, we have to take into account the possibility of dynamical mass 
generation. Generally the presence of masses is similar to hard momentum scales in that the hard momentum region 
of loop integrals can be important. The extended power counting analysis that takes into account all these aspects 
has been presented in section [ll] where the possibility of dynamical mass generation was linked to the IR exponents 
of the propagators via the auxiliary symbols in eq. (13 1 for the gluons figi and ghosts /ig/j. 



In the presence of several different scales a mere power counting of anomalous IR exponents is not sufficient anymore 
but the canonical scaling of the integrals, propagators and vertices has to be considered. In particular it is possible 
that tensor structures of a vertex involve hard momenta and do not scale with the soft momentum. Correspondingly, 
it is also necessary to discriminate between bare and dressed vertices in this context. Therefore we will first assess 
the canonical scaling of the appearing vertices in detail. Let us start with the bare ghost-gluon vertex which in our 
conventions depends only on the outgoing ghost momentum. If this momentum is soft the canonical scaling has to 
be taken into account independent of the size of the other momenta and vice versa. In addition, when the vertex is 
connected to an internal gluon propagator the tensor in the direction of the gluon momentum is canceled due to the 
transversality of the gluon propagator in Landau gauge 



q'^D^,, (p) =(^q- ^pj D^, {p) = q^^D,,, (p) . (25) 

The dressed ghost-gluon vertex has two independent tensor structures that can be chosen as arbitrary linear combi- 
nations of the 3 external momenta. If only one external momentum is soft there is a tensor structure that depends on 
hard momenta which will dominate as long as there are no cancelations. Therefore a dressed vertex has generically 
no canonical scaling whenever there are hard scales involved. However, because of the above transversality this can 
be changed when there is only one hard external scale - as is the case in propagator integrals - but not if there are 
two independent ones. 

The 3-gluon vertex has many tensor structures but it turns out to be sufficient to analyze the tree level tensor 

(ro);^; (P: 9, r) = -tgr"^ ((p - q)^ S,, + {q ^ r) ^ 8,, + (r - p)^ 8,^ . (26) 

When only one of the momenta is soft the tensor is in general of the order of the hard momenta that dominate soft 
contributions. Other possible tensor structures that depend only on soft momenta are likewise subleading compared 
to the tree level tensor and the canonical scaling is again not present. Gluon transversality cannot directly change this 
here since the above tensor structure involves the metric tensor which couples the two attached gluon propagators in 
a loop directly. However, there are cases where the Bose symmetry of the vertex can lead to additional cancelations 
as discussed below. 

Let us now apply this to the individual diagrams to obtain the corresponding equations for the IR exponents. The 
decomposition of the full ghost DSE reads 
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c\ 



Since it depends only on a single momentum the dominant loop momenta that contribute in the loop graph are 
normally of the order of the soft external momentum scale k ^ p but when masses are generated dynamically also 
much larger momenta k ~ ^qcd 3> p can contribute and could in principle be relevant due to the positive mass 
dimension of the graph. In particular, because of the possibility of kinematic singularities of the vertices one obtains 
a nontrivial contribution from large loop momenta that involve the IR exponent of the ghost-gluon vertex in the limit 
that the ghost leg becomes soft. Although the integration measure for this contribution from hard momenta does not 



scale with the soft external momentum, it is suppressed due to cancelations owing to the gluon transversality eq. ( 25 1 
which introduce a canonical scaling part for the ghost-gluon vertices. This immediately rules out the possibility of an 
IR constant ghost propagator since the leading contributions from finite modes that could produce such an IR ghost 
mass nigh (p) 3> p are canceled in the ghost integral and only subleading contributions remain. Because of the absence 



of such a mass = 1 and there is the possibility of an additional suppression of the propagator eq. (12) given by 
figi in case the gluon propagator is likewise massless. The corresponding equation for the anomalous dimension of the 
ghost reads therefore 



-Sgh + l = min (l; 6^g + 5gh + Sgi + 1, (5^^ + 1 + ^lgl) , (27) 

where here and in the following we separate terms arising from different graphs by semicolons and those arising from 

different kinematic regions of the same graph by commas. 

The gluon equation truncated to IR leading 1-loop terms is given by 



-1 IR 

s 



s 




h 




As in the ghost case there can be contributions from the hard part of the loop integral. In particular, in this case 
there is no suppression due to transversality as in the loop correction of the ghost propagator and therefore an IR 
gluon mass is not ruled out. The absence of an IR ghost mass implies that according to eq. (12 1 the contribution 



from the hard region of the ghost loop is suppressed by p^/fc^ so that this correction cannot induce an IR gluon mass. 
This is different for the gluonic corrections that can contribute once a gluon mass has been dynamically generated 
but are likewise suppressed otherwise. Since the external scale does not enter the tadpole correction at all, it is not 
necessary to distinguish between different kinematic regions in this case and as long as the there is no scale induced 
this contribution is canceled identically. Correspondingly, the gluon equation reads 



-Sgi + l= min (l; S^g + 2Sgh + l,Sfg + l; + 2Sgi + 1, ,5f^ + figi ; figi) , (28) 

The leading contributions in the ghost-gluon DSE are given by the two triangle graphs that are analogous to the 
Abelian and non-Abelian diagrams in the quark-gluon vertex. As discussed in section |ll] when scales of different 
order of magnitude are involved, the decomposition of the loop integrals yields in addition to a hard contribution 
generally two contributions from inequivalent kinematic regions of the loop integral shown in fig. [TO] where different 
internal momenta become soft. These correspond to inequivalent ways to route the large momentum through the loop. 



E.g. for hard and soft external momenta p and q, respectively, and an assignment of the loop momentum k such 
that the first IR relevant contribution in fig. 10 arises from soft loop momenta k q, the second one arises from hard 
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Figure 13: The IR leading part of the decomposed equations for the ghost-gluon vertex in the uniform (top), soft-gluon (bottom 
left) and soft-ghost limit (bottom right). 

loop momenta in a correspondingly narrow momentum interval k^p^q-^p^k. When assessing the counting of 
the vertex, each of these kinematic contributions of the initial loop integral could dominate and has to be taken into 
account separately. Using the graphical representation depicted in figs. [TO]and[TT]the ghost-gluon DSE in the different 
kinematic limits is explicitly represented in fig. [13] We should stress that there is no double counting involved here, 
since the figure only visualizes the different distinct IR-sensitive regions of the loop integral that arose from the exact 
decomposition discussed in section |ll] and appendix |B] 

From the decomposed DSE we obtain the IR counting of the individual contributions to the ghost-gluon vertex in the 
uniform limit 

+ ^ = min (^-■25^g^ + 25,^,+5,l + \,5fg+5fg+l+^,gV,5%^^ ■ (29) 

Here in the minimum function the first element is the bare vertex, the next two arise from the regions of soft and 
hard loop momenta of the "Abelian" graph whereas the final two are the contributions from the corresponding regions 
of the "non- Abelian" graph. In the contributions from hard momenta the 1 arises again from cancelations due to the 
gluon transversality. In the limit that only the gluon becomes soft we find 

5l\ = min (^)-5l + 5f^+25,,, + \,25f^+5,i + l,5fg;5-,^+5^^^^^^ . (30) 

The additional suppression of the last term from the hard momentum region of the "non- Abelian" graph arises from 
the combination of gluon transversality and the Bose symmetry of the 3-gluon vertex. Choosing as the basis for the 
dressed ghost-gluon vertex the outgoing ghost- and the gluon- momentum, Tgg = Ap^^* + B pgi, with corresponding 
dressing functions A and B, the contraction of the dressed and bare ghost-gluon vertex with the corresponding gluon 
propagators in the loop gives, respectively, 

(^A-{k + h + s)^ + B-ik + s)^) 1^5^^ _ ^^±A(|±fl^^ ^A-(^h^-'^k^^+Ois) , 

+ s)s (ds, -'^^^hp-'^kp + Ois) , (31) 

where h is the incoming hard external ghost momentum, s the soft external gluon momentum and k the loop momen- 
tum. Because of the antisymmetry of the color part, the Bose symmetry of the dressed 3-gluon vertex implies that it 
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Figure 14: The IR leading part of the decomposed equations for the 3-gluon vertex in the uniform (left) and soft-gluon limit 
(right). 

is antisymmetric with respect to simultaneous commutation of momenta and Lorentz indices. The appearance of two 
identical momenta up to O (s) corrections shows then that the leading term vanishes and the loop scales actually as 
O (s). The corresponding equation in the soft ghost limit follows similarly from the power counting analysis 

= min (^■,5l+6ll+5,,,+5,i + \,25f;+5,>, + 2,5f^-5l + ^^^^^ . (32) 

Finally in the considered truncation there is the equation for the 3-gluon vertex where the contributions are given 
by the ghost and gluon triangles. The equations for the two different kinematic limits are shown in fig. |14[ In the 
uniform limit the power counting yields the equation 




where the 1/2 in the contribution from hard modes of the gluon triangle arises since the leading term depending 
only on the loop momentum is odd and vanishes in the symmetric integration. The soft-gluon limit requires some 
more care. Here, the IR exponent of the ghost triangle seems to depend on different factors like the loop routing, the 
definition of the bare ghost vertex or which of the three vertices is taken bare in the DSE. All these factors seemingly 
determine whether the appearing bare ghost-gluon vertex scales canonically. To see that this is actually the case 
independent of all these conventions, it is important to remember that the 3-gluon vertex is totally symmetric. This 
means that all three ways of assigning the external momenta and indices to the external legs of the graph yield the 
same result. With the standard convention that the bare vertex is proportional to the outgoing ghost momentum, 
this momentum is soft for at least one of the three configurations. Since the IR exponent obtained by the power 
counting analysis can overestimate the degree of divergence of a given graph when there are cancelations but it cannot 
underestimate it, the bare ghost-gluon vertex features indeed a canonical scaling. Correspondingly one obtains here 



= min (^■5^^^+5f^+25,H + l,5f^+5f,+5,h + \,5f^^^ ■ (34) 

The fourth element in the minimum function gives a constraint on the two vertices in the soft gluon limit 6^g < S^g. It 
is straightforward to see that the reduced system of the three equations for the soft divergences alone is consistent with 
the uniform solution discussed in the last section if the corresponding IR exponents are taken as given. Though, here 
we want to discuss the general case where the additional kinematic divergences couple back and could thereby change 
the uniform solution. To this end we study the full, coupled system of seven equations obtained by the consideration 
of possible kinematic divergences. 

We want to stress that, analogous to the case discussed in detail in section [Till the IR exponents on the left-hand side 



in the vertex eqs. ( 29 1 , ( 30 ) , ( 32 1 and ( 34 1 appear linearly also on the right-hand side. The corresponding inequalities, 
e.g. Sgg + 1/2 < dgg+S!jg+Sgfi + 25gi + lj2 from eq. (29), can be fulfilled independently of the values of the IR exponents 
on the left-hand side, since they drop out. Yet, they provide constraints on the other IR exponents. From a purely 
mathematical point of view these DSEs could be trivially fulfilled for a whole range of values for the corresponding 
IR exponent if these linear terms would dominate. Yet, such linear terms provide no non-linear feedback and cannot 
alter the degree of divergence compared to the tree level term in the DSE that is not IR enhanced. Correspondingly 
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from a physical point of view these terms are irrelevant for the IR solution. 

As in the uniform case the system of IR equations is again constrained by necessary conditions for a stable skeleton 
expansion. In addition to those in the uniform limit eqs. (22 1 there are corresponding constraints from the graphs in 



fig. [12] when the loop momentum is in the IR regime but one of the two connected propagators has a hard momentum 



or when the momenta of both connected propagators are hard 



(35) 



2'5fg + <5„, + 1 > 



5gl + l>0, 

^gi + 1 > , 
'gi 



(36) 



In principle there may also be constraints from the hard part of the added loop, but in this case the extension can 
chan ge t he counting fax away from the insertion and therefore they do not ta ke a simple form. Using the constraints 
eqs. (22 1 , (35 1 and (36 1 as well as 5fg < S^g, the system of inequalities, eqs. (27)-(34l, reduces to: 



-6gh + 1 = min (1, (5;; + Sgh + Sgl + 1, S^g^ + 1 + flgl) , 

-5gl + 1 = min (l, S'^g + 25gh + 1, ^3",, + 28gl + 1, Sfg + flgl^ 

^n^(^l,Sll+Sfg+l + t,gl^ , 



"gg 

ggi 
"gg 



mm I l,2<5g"g+3<5,, + i,25f; + i 
iTvm{Q,25fg+5gh + 2) , 
mm 



in(^^5-gg+5f:g+25gn + \,5fg+5fg+5gn+2,5fg + , 
min {Q,5-gg+5fg+25gu + l,6fg+6fg+5gn + \,5f)i . 



(37) 



In the uniform case the non-renormalization of the ghost-gluon vertex was the cornerstone that allowed one to solve 
the corresponding system. We therefore start with this vertex in the soft ghost limit described by the equation for 
Let us assume for the moment that the vertex is singular so that the second term dominates. In this case we 
would have = —5gh ~ 2, but this would yield a direct contradiction via the last term of the ghost propagator 
equation. Therefore we have instead (5|g = and obtain additionally the weak constraint 5gh > —2. In contrast to 
the conformal case the remaining system allows two qualitatively different solutions: 

Decoupling solution: With any renormalization prescription that does not change the above system for the IR 
exponents, as in the conformal case discussed above the equation for the ghost propagator as it stands has only the 
trivial solution 5gh = 0. To find the solution of the residual system in the present case, we insert the remaining 
possible solutions for the 3-gluon vertex in the soft gluon limit into the other equations which yields 



-&gi + l = min (1, ^,gl) , 5-g = 6fg = 5^ = = . (38) 

In addition to the trivial fixed point, there is here also the possibility of a massive gluon Sgi ^ 1. Therefore, in contrast 
to the conformal analysis where it was assumed that all integrals are dominated by scales of the order of the external 
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momenta, hard modes of the order of the induced scale Aqcd can dominate the loop integrals in the equation for the 
gluon propagator and induce an IR gluon mass. This alternative scenario is reminiscent of the early work [24_ and has 
recently also been suggested in |21fB5] on the level of the propagator equations, whereas the present analysis shows 
that it is also consistent with the vertex equations. This solution is also found if the ghost propagator is renormalized 
to a finite value at vanishing momentum [23 . 

Scaling solution: Alternatively, it is again possible to use the renormalization introduced in |1H IT3] as discussed 
in the last section which removes the tree-level term in the ghost equation. Because of the scale independence of the 
ghost-gluon vertex in the limit that a ghost momentum vanishes the last term in the ghost equation is either directly 
subleading in case the gluon propagator remains massless or can otherwise be removed in the renormalization process 
together with the tree-level term. Thereby the ghost equation becomes unique 

Sgh = -l{Ssl+S^,) (39) 
and inserting this expression for the ghost exponent, the gluon equation becomes trivially fulfilled 



-5gi + l = min (l, -6gi + 1, + 25gi + 1, ,5f^ + A^g;) , (40) 

which shows that the ghost-loop has to be the IR leading contribution. It does not provide any constraint on Sgi, 
leaving the free parameter k = dgi/2 > in the solution as is known already from the uniform case. We note that the 
triviality of this equation is qualitatively different from what we found in the case of the equations for the kinematic 
divergences above, since the equation is non-linear and can therefore indeed self-consistently generate a non-trivial 
solution. The equation provides two additional weak constraints on the vertices > —6k, S^g > 1— 2k— /ig/. Inserting 
the above expression for the ghost exponent in the first of the constraints from the skeleton expansion in the uniform 



limit eqs. (22 1 yields S^g > which makes the corresponding equation for the exponent of the uniform ghost-gluon 
vertex unique, so that 5gg — 0. The residual system reads now 

Ssg + 1= min (^-3k + i, 2,5f^ + ^ + t^gi 

5fg = min (^0,-2K+|<5f;-K + 2,<5f; + \ 
Sll = min f 0, -2^+1, Sfg-K+^,S A . (41) 



Using the constraint (Jg^ < 1 — 2k from the last of these equations in the third of eqs. (35 1 yields the constraint 
K < 3/2, so that a solution can only exist in a bounded region of the IR scaling parameter k. The system can then be 
solved by inserting the possible solutions of the ghost-gluon equation in the soft-gluon limit in the equation for the 
3-gluon vertex in the same limit. Under consideration of this bound the equation reduces to S^g — min (0, — 2k-|-1). 
Correspondingly, the scaling of the 3-gluon vertex in the soft-gluon limit depends on the value of the IR parameter k. 
For K > 1/2 it is singular Jf^ — 1 — 2k whereas for k < 1/2 it is not S^g = 0. Using these results the remaining equation 
becomes trivial and the scaling solution of the system is obtained. Similar to the 3-gluon vertex the exponent in the 
soft-gluon limit of the ghost-gluon vertex min (O, | — 2k) also depends on k, but here a divergence would arise only 
for K > 3/4. Finally, the scaling of the 4-gluon vertex in the uniform limit that satisfies the DSE system is obtained 
from the corresponding DSE fig. |6] as before and yields the known result S^g = —4k whereas the determination of 
the corresponding kinematic divergences requires a more detailed study. In summary, the IR fixed points for the two 
qualitative distinct solutions are given by Table [Tj 

These solutions fulfill all constraints that appeared in the course of the evaluation and present therefore refined IR 
fixed points of Landau-gauge Yang-Mills theory. Several remarks are in order at this point: 

• In this work we restricted our analysis to the IR behavior of Yang-Mills theory. The physical relevance of such a 
study stems certainly from the fact that Yang-Mills theory presents the gauge sector of QCD and might therefore 
provide important insight into qualitative properties of the strong interaction. Whereas an independent study 
of the gauge dynamics is by definition sufficient in the quenched limit [B] , in the case of dynamical QCD it is 
not a priori clear that the quark dynamics does not affect the IR fixed point structure of gauge Green functions 
obtained here. Yet, a corresponding recent IR analysis of QCD [7J, based on the general methods developed in 
this work, shows that the gauge sector is totally unaltered by the quark dynamics and thereby strongly supports 
the relevance of the results presented here and in other studies of Yang-Mills theory. 
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Table I: The IR exponents for the leading Green functions of the IR fixed points of Landau gauge Yang-Mills theory within the 
two possible IR scenarios. The scaling analysis restricts the parameter to be only positive k > 0, and yields rather weak upper 
bounds. (Note that these bounds on k are weaker than the one that had been erroneously given in a preprint version of this 
article. Stronger bounds 0.5 < k <1 are, however, supported by analyses of the actual loop integrals ^12).) 



• In contrast to the case of the conformal analysis discussed in the previous section, where the chosen boundary 
condition merely excluded other solutions, the decoupling solution obtained for a generic condition is qualita- 
tively different from the scaling solution obtained for a choice that had an unbroken global BRS charge [23J. 
Whereas in the scaling solution the ghosts are strongly IR enhanced resulting in divergent gluonic vertices, in 
the decoupling solution neither the ghosts nor the vertices are anomalously enhanced. This strongly suppresses 
any IR dynamics mediated by the gluons in the ratio p = /rn?gi where nigi is the finite IR limit of the gluon 
polarization. Thereby it is easy to see that in both cases the leading contribution to a Green function is given 
by the ghost dynamics. In the decoupling case the IR exponents in Table |l] show directly that the leading term 
in the skeleton expansion of a general vertex with n ghost-pairs is anomalously suppressed by whereas purely 
gluonic Green functions scale canonically. In addition to the given decoupling solution there might be further 
IR fixed points where even the vertices decouple and become IR constant [33J. 

• Recent lattice simulations on large lattices ||3j show a gluon propagator that does not feature a decrease in the IR 
and a ghost propagator that is basically not IR enhanced and thereby favor the decoupling scenario. Moreover, 
it has been argued that this is probably neither a finite volume [34J nor a statistical effect [35J. However, there 
seem to be issues with Gribov copies ||36j |37j and discretization effects [38j in these analyses. Such effects could 
shadow a potential scaling behavior in the deep IR, so that the scaling scenario cannot be excluded with the 
present data. Another view is presented in f39', where it is suggested that the different boundary conditions 
|23| correspond to distinct residual gauge fixing conditions. In principle Gribov copies could be an issue in the 
DSE system as well and the constraint to the fundamental modular region that is free of Gribov copies could 
change these equations and their solution structure. However, it has recently been shown explicitly [40 that at 
least the restriction to the first Gribov region [5^ using the Gribov-Zwanziger action [H] does not affect the IR 
fixed point structure of the scaling solution. 

• The IR exponents given in Table |I] determine only the anomalous scaling laws for the most singular tensor 
parts. The scaling of the full Green functions involves also the canonical scaling dimension incorporated in the 
tensors. In particular it is possible that some dressing functions are more IR singular than the leading dressing 
function, but their contribution to the vertex is nevertheless subleading since it is additionally suppressed by 
the canonical scaling of their tensor. As we show in [5(1 this is indeed the case for the ghost-gluon vertex 
which features more structure than the above result suggests. Instead, a soft-gluon singularity appears in the 
form factor of the longitudinal tensor that is additionally suppressed by the gluon momentum in the tensor and 
actually IR vanishing, whereas the tree-level tensor is entirely IR finite and presents the IR leading structure. In 
order to reveal such subtleties in our power counting analysis we would have had to include different anomalous 
dimensions for the different tensor structures. Since we present an explicit analytic solution for the IR limit 
of the 3-point vertices in [20) we refrained here from such complications. In [20j we show explicitly that the 
transverse part of the 3-gluon vertex does not have leading tensors. The analytic result obtained there supports 
this observation. It also explains why the kinematic singularity for the three-point vertices obtained in the 
second work of ref. [TB], where only the transverse part was considered, is lower as the one given in tab. [ij 

• The kinematic singularities do not alter but merely extend the previously know uniform scaling fixed point. Yet, 
the soft singularities restrict the range of possible k- values from the mere positivity requirement in the conformal 
case to the bounded interval < k < 3/2. The trivial solution is included in the scaling solution for k = 0. It 
is (up to logarithmic corrections) realized in the UV regime of the theory characterized by asymptotic freedom 
and it is clear from the perturbative /3-functions that this solution can hardly be a stable IR fixed point, too. 
The best currently known value for the IR scaling parameter is k sa 0.5953 [121 IT^ obtained from an analytic 
IR solution of the integrals in the DSEs for the propagators. 

• Interestingly below k = 1/2 the kinematic singularities entirely disappear. The latter value is a special case 
since the gluons show effectively a massive behavior but the ghost is in contrast to the corresponding decoupling 
solution still strongly divergent. For k > 1/2 the gluon propagator should vanish with the small exponent 2k — 1 
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which is precisely the negative of the exponent for the mild kinematic singularities of the 3-gluon vertex found 
here. In four dimensions such a small exponent naturally poses a huge numerical challenge and is not observed in 
current studies |3] , but in lower dimensions a corresponding decrease has been clearly confirmed [H] . Similarly, 
it is not surprising that the predicted kinematic singularities have not been seen so far in present vertex studies 
|43| which are numerically even much more challenging than those for the propagators. In j20) we proposed a 
suitable tensor contraction for the three-gluon vertex to investigate the existence of kinematic singularities on 
the lattice. 

It is crucial that the ghost-gluon vertex is finite when only a ghost momentum vanishes. This result follows 
immediately from the corresponding "un-decomposed" DSE which contains only a single loop graph involving the 
connected (instead of IPI) ghost-gluon scattering kernel [H] . By transversality this graph is directly proportional 
to the external momentum and leaves only the tree level part in the IR limit in accordance with the non- 
renormalization of this vertex. We point out that there is, however, no corresponding argument when the gluon 
momentum vanishes. 

The obtained divergence when only a single gluon momentum vanishes naively seems to be problematic for 
several reasons: First of all it seems to induce an even stronger singularity in the ghost-gluon vertex in the 
uniform limit from hard loop momenta. As pointed out above though in this case the transversality in Landau 
gauge prevents this and instead makes this contribution strongly subleading. This is also in accordance with the 
two different versions for the ghost-gluon vertex. Since a dressed 3-gluon vertex is present only in the first one, 
the two versions would be inconsistent if the kinematic-divergence of the 3-gluon vertex would alter the degree of 
divergence of the full vertex. Secondly, naively there seems to be a huge problem with the soft gluon singularity 
in the 3-gluon vertex. First of all it arises directly from the ghost loop integral with dressed propagators. But 
once induced, it seems to arise in addition also in dressed vertices whenever the external momentum becomes 
soft and totally independent of the loop integral. This would enhance the divergence in each iteration and make 
it more and more divergent. As seen explicitly in the above analysis the reason why this is not the case is 
that the hard region is additionally suppressed and that both of these different singularities arise from distinct 
regions of the loop integration and thereby cannot amplify each other. 

As found from the analysis above, all Green functions in the scaling solution include a graph that does not 
involve singular vertices. In particular, the IR dominant ghost loop correction to the 3-gluon vertex induces 
the soft-gluon divergence entirely due to the enhancement of the ghost propagator so that the appearance of 
kinematic divergences is a direct consequence of the ghost dominance of the uniform solution. This allows 
one to capture the qualitative IR behavior of the vertices in a semi-perturbative scheme that involves dressed 
propagators but employs bare vertices. This approximation is used in a companion article and allows a complete 
analytic solution in terms of hypergeometric functions L2Q_ . Although the analysis of the kinematic divergences of 
the 4-gluon vertex is more complicated and requires a detailed analysis, we note here that the semi-perturbative 
contribution from the ghost loop yields Jf^ = 1 — 2k and = 1/2 — 3k which due to the ghost dominance is 
expected to be the leading contribution. Incidentally this is exactly the scaling of the 3-gluon vertex in the 
soft-gluon and uniform limit whereas the uniform limit of the 4-gluon vertex 154^ = —4k is more divergent by 
yet another — 1/2— k. These results suggest the conjecture that the full scaling of a gluonic correlation function 
with rig independent soft external momenta should be r„^ ~ (P^)^ {n^+i){i/2+K) ^ 

Inserting the results for the scaling solution in Table [T] in the constraints for the skeleton expansion involving 
hard momenta we find that the extensions count as 3/2 — k and 3 — 2k when there are one respectively two hard 
propagators in the extended graph. With the above limits for k this shows that these extensions are strongly 
suppressed and correspondingly such extensions do not have to be taken into account in the skeleton expansion. 



In contrast, inserting the results in the constraints for the skeleton expansion in the uniform limit eqs. (22 1, it 



is clear that all of them are saturated and correspondingly all orders in the expansion scale equally \}Aj . Note at 



this point, that we did not assume by our constraints eqs. (22 1 that the skeleton expansion strictly converges, 
but only that it is not explicitly divergent. In general, the skeleton expansion could be an asymptotic series as 
suggested by the IR scaling. Therefore, the whole tower of such graphs had to be resummed which could in 
principle change the IR scaling. For instance it is well known from standard resummed perturbation theory that 
the resummation of perturbative logarithms yields a power law scaling with an anomalous exponent. However, 
even such logarithmic divergences can be invariant under resummation, as e.g. found for the non-Fermi liquid 
corrections in dense QCD |33]. In the current case a resummation seems to be impossible in full generality, 
anyhow, but the decisive difference is that the graphs that are resummed already feature power law scaling 
and therefore we expect that they are indeed invariant under resummation. Furthermore, it has been shown in 
|15l [?fl H5] that the uniform IR solution can be obtained independently of the skeleton expansion and the same 
will be shown for the more general case as discussed in a forthcoming publication [TJ . 
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• More generally the validity of the skeleton expansion is closely linked to the existence of any finite truncation 
of the DSEs and eventually to the concept of locality. To see this, note that an explicitly divergent skeleton 
expansion suggests that there is no finite approximation to describe higher order Green functions by lower ones 
or in particular only in terms of the primitively divergent correlators. This would mean that the higher Green 
functions include important physics that is not yet included in the lower Green functions and that is required 
to properly describe the system. In particular these higher Green functions in turn significantly influence the 
solution of the leading Green functions within an explicit analysis of the dynamics in this case. Any truncation of 
the system with a local effective action that includes only a finite number of terms would thereby miss the main 
physics and instead a non-local effective action with an infinite number of terms in the local fields is necessary. 
Such a non-local situation means that the description in terms of local degrees of freedom breaks down and the 
system can probably be more conveniently described by non-local degrees of freedom with a finite number of 
terms in the action. Now one could argue that, since one is only interested in the behavior of the leading Green 
functions, as far as the appearing higher order Green functions in the corresponding DSEs are known (e.g. the 
5- & 6-point functions in Yang-Mills theory) the system for the primitively divergent Green functions may be 
solved independent of these problems. The important point to realize is that in such a case there is no way 
to obtain them, neither in the context of functional techniques like DSEs nor in any other scheme like lattice 
gauge theory since even there the number of lattice points corresponds to the highest Green function that can 
be realized and that could contribute to the dynamics of the system in this truncation. Correspondingly, in 
contrast to the usual situation where one expects that the precise form of the unknown higher Green functions 
becomes irrelevant at sufficiently high order and thereby a sensible ansatz should be given by general arguments 
like symmetry restrictions, this is intrinsically not fulfilled in case of an explicitly divergent skeleton expansion. 
Finding a reasonable ansatz for these higher order Green functions is thereby equivalent to guessing the correct 
solution for the primitively divergent Green functions in the first place. Finally, even if we would simply by 
chance guess the precise solution for the required higher Green functions to solve the system for the lowest 
correlators like the propagators, the above line of reasoning shows that these local Green functions would have 
nothing to do with the actual physics of the system since these degrees of freedom should not be suitable to 
describe the system in case of a divergent skeleton expansion. From this point of view we regard the existence 
of a stable skeleton expansion as a rather physical requirement for any analysis in terms of underlying local 
degrees of freedom. 



VI. CONCLUSIONS 



We have presented a power counting formalism for the analysis of the non-perturbative IR behavior of field theories 
in the general case when finite scales are present. This allows the inclusion of kinematic divergences and dynamical 
mass generation in the power counting analysis. In this framework we studied the IR regime of Landau gauge Yang- 
Mills theory in more detail and found that the fixed point structure is more diverse than previously assumed. As a 
general result that does not rely on any approximations we find that the DSEs directly exclude a class of IR fixed 
points where the IR strength arises directly from the gluon dynamics. Instead there are two other qualitatively 
different solutions for the fixed point behavior. 

In the scaling solution the IR regime is strongly dominated by the ghost dynamics, as predicted e.g. in |1H416|, |47]. The 
structure of the scaling fixed point established in these studies, however, has to be amended by additional kinematic 
singularities. The presence of these singularities is not only consistent with the uniform scaling rules but has both 
conceptual and quantitative impact on the structure of this solution. In a companion article [2Q^ we present detailed 
analytic results for the 3-point vertices that give the complete kinematic dependence and show precisely the same 
kinematic divergences found here by pure power counting arguments. The knowledge about the leading dynamical 
contributions obtained here in combination with the analytic results for the 3-point functions should allow to give an 
improved value for the IR exponent k in an approximation that treats the 3-point vertices dynamically. This should 
include the main dynamical contributions for a precise prediction of this pertinent parameter. 

The non-perturbative analysis of kinematic divergences established here is even more important in the case of QCD. 
There non-perturbatively enhanced, strong kinematic divergences of the quark-gluon vertex in the corresponding 
scaling solution can provide a description of crucial aspects of QCD, like the linear rising potential between static 
color sources in the quenched theory [6J and its long-range screening in the dynamical case[7| as well as the U(l) 
anomaly f48]. 

Depending on the renormalization prescription there is another decoupling fixed point were the gluon acquires a mass 
and decouples |21ff55] . Current lattice results for the propagators clearly favor this possibility Yet, we find that, 
aside from possible inherent problems of these analyses, cf. |36H38j . in the decoupling scenario the gluonic vertices 
remain bare in the IR limit and thereby this fixed point shows no infrared enhancement at all. Correspondingly it does 
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Figure 15: The generating ghost-DSE. Thin and thick lines and sohd dots represent bare and proper propagators and vertices. 
Open dots represent explicit fields and the black square separating the two different lines denotes the corresponding off-diagonal 
component of the super-propagator. 



not provide a description of the above vital aspects of QCD in terms of the Green functions of the fundamental local 
degrees of freedom alone. Despite its apparently simple form it would instead require a more complicated description 
of the QCD vacuum in terms of other degrees of freedom. 
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Appendix A: Graphical derivation of Dyson-Schwinger equations 

The Dyson-Schwinger equations for the vertices given in the main text were derived via an algorithmic method 
presented in [26 . In this appendix we sketch this method to derive Dyson-Schwinger equations for general correlation 
functions from the corresponding equation for the 1-point function which we will in the following refer to as the 
"generating DSE". They involve except for one bare vertex only proper correlation functions and are given for the 



ghost and gluon in figs. 15 and 16 respectively. As shown in detail in [26' within a convenient superfield formalism all 
other DSEs can be computed from these equations via the replacement rules in fig. [17] where the double lines stand 
for superfields that include all elementary fields in the theory. As the final step all expressions have to be evaluated 
at their vacuum expectation value which corresponds graphically to replacing all super-propagators and vertices by 
the irreducible ones in all possible ways that involve only physical propagators and vertices in accordance with the 
symmetry of the action. This generally removes many graphs, so that when performing this procedure it is useful to 
take into account beforehand to what order the DSEs shall be computed to neglect any unphysical terms that would 
already vanish anyhow during the extension steps. 

In the special case of non-Abelian gauge theory the super-multiplet containing gluon and (anti-)ghost fields is given by 



(j) = {A, c, c) and denoted by curly and dashed lines respectively. Using the general graphical replacement rules fig. 17 
yields in a straightforward way the equations for the leading correlation functions figs. ^to[6] The proper symmetry 
factors arise simply from different ways of obtaining the same graph in the above replacements steps. The arising 
super-propagators in the 2-loop term of the generating gluon-DSE fig. [16] only become relevant for the equation for 
the ghost-gluon vertex and corresponding higher order correlation functions. For the derivation of other correlation 
functions they may be replaced by ordinary gluon propagators. It is important to note, however, that it is not possible 
to derive the equation for the gluonic vertices by additional derivatives of the propagator equation fig. [2] since one 
would miss additional loop contributions from the next to last graph in fig. [16] that guarantee the symmetry of the 
corresponding equations in a perturbative approximation. Finally we remark, that since there are no fundamental 
quark-ghost interactions in QCD the graphical expressions for the generating quark DSE and the equations for the 
leading correlation functions are identical to the corresponding ghost equations displayed in fig. [15] respectively figs. 
[T][3]and[4] 
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Figure 16: The generating gluon-DSE. The Ihs in the generating equations represents the field derivative of the effective action. 



Figure 17: The replacement rules for the generation of DSEs for general correlation functions. One replacement step corresponds 
to a functional derivative that generates the DSE for the corresponding correlation function with one more external leg. In 
such a step one of the objects on the Ihs is replaced by the corresponding expression on the rhs. Namely, an explicit field is 
simply removed. A general off-shell propagator is extended introducing a new proper 3-point vertex where each double line can 
stand for any field in the super- multiplet. A general proper vertex represented by the thick dot (which can already have any 
number of external legs) is simply extended via attaching another leg. 

Appendix B: Decomposition of the loop integrals 
1. Massless integrals 

In this appendix we demonstrate explicitly how the loop integrals in the DSEs can be decomposed into different 
parts when both soft and hard external scales are present as the IR limit is taken. A general 1-loop graph is shown 
in fig. 18 The external momenta qn are assumed to be divided into two disjunct subsets {qn} = {si} U {hj}. The 
Si are soft and the hj hard, so that in the IR limit \si\ <^ \hj\ Similarly the linear combinations that can be 

formed out of these momenta and arise in the integral lines of a general loop graph, 



(Bl) 



with Q G {—1,0, 1}, can likewise be divided into corresponding subsets {Qn} = {Si} U {Hj} with \Si\ <ti \Hj\ 

The integral I for a particular term in the general loop integral Q, eq. (|6|, arising from the expansion of the product 

of decomposed vertices eq. (p|) reads then 



1 = 



(2^ 



lC{i.} ({(fc + Qf}, {Q'}) n (fc + Qnf'"-'^ UpI ({(^ + Q)'} ' {Q'} 



(B2) 



where n runs over the propagators and a over the vertices in the graph. Here we skip the indices when denoting the 
corresponding sets arising as arguments of the functions K. and p defined in eqs. ([6| and ([2|. The contributions to the 
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Figure 18: Exemplary momentum routing for a general n-point function. Here the arrows denote the momentum flow and the 
particular linear combinations arising in the diag ram are given by Si — 0, Hi — /ii, S2 — hi — /12 hi, /i2; Ss — hi — h2 si, 

vertex functions with different IR limits in the decomposition eq. (|2| can be divided into a subset where the power 
law contains hard external momenta and another that does not 

J (27r) J ^ . ^ 

■ UpI ({(fc + Sf}Y' UpI ({(^ + Sf},{{k + Hf})'' . (B3) 

a b 

In the IR limit one can choose some intermediate scale &, with S <^ Q <^ H where S = max {\Si\) and H = min (|i?j|), 
which can e.g. be chosen as = V SH. Then the loop integral can be split into two parts 



1 = 



' pdk + rdk] fc-i / (• • • ) n + s^f '^"' IK^ + H,f ''"' UpI (• • • UpI (•••)' 

v-'o Je J J (27r) ^ ^ a b 



= !<+!>. (B4) 
Consider first X<, where the factors involving the Hj ^ can be expanded in k, 



i<=ic,,.AW)UHr'-''Upiimf 

j b 

Here the analyticity of the pi as well as the partial analyticity of K, has been used. Note that possible non-analyticities 
of /C{;^} at hard momenta are no problem as far as the external momenta qi of the Green function themselves do 
not represent a singular configuration at hard external momenta. As discussed in the main text such a physical pole 
occurs in Euclidean space away from the real axis and directly at such a point the solution of the Dyson-Schwinger 
equations is not well defined anyhow and we can therefore simply exclude this special kinematic case here. 
The integral can now be extended over all scales. The leading part is 
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(B6) 



where here and in the following the symbol O means that subleading terms are of the given relative order and 
correspondingly suppressed compared to the leading terms that are given explicitly. In the correction term the other 
factors can be expanded 



, b Je J (27r) ^ V 



J b 

Analogously one can expand the upper part of the integral 



(B7) 



■Upl(>^Y'Upi{k^{ik + Hf}f (^l + Y.0{^§)^ (B8) 

and subsequently the correction term so that their sum gives 

AI = AI< + AI> (B9) 

This integral involves no scale and thereby simply represents a, possibly infinite, constant multiplying its prefactor. 
This prefactor is in particular completely independent of soft scales and AT does thereby not affect the IR power 
laws of the corresponding Green functions. In the case of primitively divergent Green functions considered in this 
work such a constant contribution is already present in the DSEs in the form of the tree level term. Moreover, for 
primitively divergent Green functions the integrals /SX can be explicitly UV divergent. Such a divergence cancels by 



construction precisely the additional UV divergence introduced in the additional scale dependent integral via eq. ( B6 1 . 
Correspondingly, this constant contribution could be formally absorbed into the multiplicative renormalization factor 
of the tree level term, which thereby has the form to cancel all divergences in both scale dependent integrals. We 
note that e.g. in dimensional renormalization the above scale independent integral vanishes by definition. Thereby, 
after proper renormalization the initial integral is split into two finite contributions that are in the IR limit S H 
independent of the scale 9 

j b l^"'-' i a 

\ ^) j b a 

The first integral depends now only on soft external momenta, whereas the second one depends only on hard momenta. 
These integrals are generally dominated by the poles of the integrand and have to scale with the respective external 
scales. Correspondingly once the contributions from scales of different order of magnitude are separated their scale 
dependence can be analyzed by mere power counting. 

Yet, the hard integral can still be subject to soft singularities when differences of the hard momenta become soft. 
Then the formally absent soft momenta re-enter via the condition of momentum conservation. To see this, consider 
a 3-point integral with one soft and two hard external momenta. For better readability we skip in the following the 
vertex dressing functions and the kernel /C but note that the derivation proceeds analogously with these complications. 



The decomposition eq. (BIO) reads then 



j(3) ^ ^2i^+P~2) f^^k- sf''-'^ + I ^ {k + h- sf"-'^ [k + hf ^-^'^ fc2(^-l) + O f M V (Bll) 

J {2Trf J {2nf \\h\J 
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The momentum combination k + h — s in the hard part of the loop integral is again sensitive to the soft momentum 
s in a small region around the loop momentum k = —h. To conveniently study this region it is useful to shift the 
integration variable via the substitution k' = k + h in order to divide the integral at some scale G 



(B12) 

In complete analogy to the general case above (using si — s, S2 — 0, hi — ~h and S.^^i — a, Si^2 — P- — 7) this 
integral can then be decomposed once more so that the full 3-point integral is effectively divided into three distinct 
integrals extending over all scales 



ji3) ^ f^2(a+0-2) (fc_,)2(,-l)^^2(,-l) _ ,)2(„- 1) ^2(^-1)^ [^0^+0-2) ^k + hf^-'^+O ( . 

J {2Trf J {2Trf J{2^f \\h\ ) 

(B13) 

The first two integrals are, up to O {\s\ / \h\) corrections, completely IR dominated whereas the third one is dominated 
by hard scales. This division is graphically displayed in fig. [TOjin the main text by distinct loop graphs with color- 
coded soft and hard propagators. The case of 2-loop contributions in the DSEs can be handled accordingly, with the 
complication that both loop integrals have to be split and the number of distinct integrals increases. In general the 
number of possible divisions depends on the n-point function, the number of loops and the numbers of soft and hard 
external momenta. After a full decomposition of the integrals there are - aside from purely hard integrals that yield 
a scale independent constant - only integrals left that depend on a given set of soft external momenta {s^}. These 
momenta can be parameterized in the form Si = pi s with a single scaling variable s and fixed, bounded momentum 
ratios pi = Si/s, \pi\ < 1 that could likewise be replaced by angular variables. In such a case the qualitative IR scaling 
of the corresponding integral is independent of the pi and can be determined by a pure power counting analysis. The 
prefactor of the IR scaling law however involves the ratios pi and requires an explicit solution of the integral |20| . 
For the cases we study in this work, the integrals that need to be decomposed depend even on a single soft scale and 
therefore trivially scale as some power of it. 



2. Massive integrals 



In this subsection we extend the previous decomposition of loop integrals to the case when mass scales are involved. 
As already discussed in the main text, in the case that all propagators in an integral are massive one can in the IR 
limit p <^ m simply expand the whole integral, e.g. 

d'k \ — i-^ « / d'^fc ^ ^ + O ("M") (B14) 

to convince oneself that the integral does not depend on the soft external scale to leading order. Even though there 
is no pole here the integral is dominated by modes of the order of the mass: the contribution from low modes is 
suppressed by phase space and the renormalized contribution from large modes is suppressed by the denominators. 
When there are both massive and massless propagators present the integral receives contributions from both soft and 
hard scales. To see this consider the massive 2-point integral of the form 

Ira= I d'^k j-^ ,, ^ . (B15) 

J {k+pfk^+m^ ^ ' 

In the limit p <^ m one can again introduce an arbitrary intermediate scale Q and divide the integrals 
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™2 7o 7 (^k+pf^K \rnJJ Jq J k^»k^ + m^\ \\k\ J J 

= \ rd'k ^+ r d'k^TY^-\ r d'k-^+oi-^ (B16) 

m? (k+p) Jo k^"k^+m^ Jo k^'^ \m J ^ ^ 

where the scale independent integral again cancels momentum independent divergences introduced in the decompo- 
sition. The first integral is dominated by soft scales and the second one by hard scales of the order of the mass and 
correspondingly the scale separation in the IR limit allows again the use of a power counting analysis. An analogous 
decomposition is possible for massive vertex integrals as well as for integrals that involve both masses and hard exter- 
nal momentum scales. In all these cases the presented decomposition results in a complete separation of the various 
external scales in different loop integrals and allows thereby the application of power counting methods. Finally we 
note that the presented method for decomposing massive integrals is likewise suitable for multi-loop corrections, like 
the 2-loop graphs arising in certain DSEs, via successive decompositions. 
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